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ABSTRACT 

I calculate the action of a satellite, infalling through dynamical friction, on a coplanar gaseous 
disk of finite radial extent. The disk tides, raised by the infalling satellite, couple the satellite and 
disk. Dynamical friction acting on the satellite then shrinks the radius of the coupled satellite-disk 
system. Thus, the gas is "shepherded" to smaller radii. In addition, gas shepherding produces a large 
surface density enhancement at the disk edge. If the disk edge then becomes gravitationally unstable 
and fragments, it may give rise to enhanced star formation. On the other hand, if the satellite is 
sufficiently massive and dense, the gas may be transported from ^ 100 pc to inside of a 10 to 10s of 
parsecs before completely fragmenting into stars. I argue that gas shepherding may drive the fueling 
of active galaxies and central starbursts and I compare this scenario to competing scenarios. I argue 
that sufficiently large and dense super star clusters (acting as the shepherding satellites) can shepherd 
a gas disk down to ten to tens of parsecs. Inside of ten to tens of parsecs, another mechanism may 
operate, i.e., cloud-cloud collisions or a marginally (gravitationally) stable disk, that drives the gas 
< 1 pc, where it can be viscously accreted, feeding a central engine. 



Subject headings: galaxies: nuclei 
disks 



1. INTRODUCTION 
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The tidal interaction between small satellites and the 
gas or particle disks in which they are embedded is a sub- 
ject of wide study in the planetary community. These 
satellites are tidally coupled via disk tides, i.e., satellites 
excite spiral density waves at the Lindblad resonances in 
the disk (Goldreich & Tremaine 1978, 1980; Artymow- 
icz 1993). As a result of these tidal interactions, gaps 
can be opened up in the embedded disks as in the case 
of the shepherding moons of Saturn's rings (Goldreich 
& Tremaine 1978) or in type-II migration (Ward 1997). 
Alternatively, if a gap does not open up, satellite may 
migrate inward rapidly due to tidal torque imbalances, 
i.e., type-I migration (Ward 1997). 

The wide applicability of satellite-disk interactions in 
the planetary community raises an interesting question 
of whether the same physics may be applicable at larger 
scales, i.e., galactic scales. In this paper, I will address 
this question by considering a very simple model, which 
illustrate the modification of the physics of satellite-disk 
interaction when applied on a galactic scale. Most no- 
tably, the critical difference is that in addition to tidal 
torques, the satellite will experience dynamical friction. 
The inclusion of dynamical friction produces a non-trivial 
effect. Namely, dynamical friction on the satellite pro- 
vides a sink of angular momentum in the system. As a 
result, the coupled satellite-disk system will continually 
lose angular momentum and sink toward the center. 

To study this physics, I consider a very simple model. 
In my model, a satellite starts out in a circular orbit at 
a large radius and sinks toward the central mass concen- 
tration because of dynamical friction on the background 
stars. Along its infall, it encounters a coplanar gaseous 
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disk or ring, which initially has a finite radial extent, r^^. 
Tides begin to couple the satellite with the disk. Because 
the satellite-disk system continues to suffer an ongoing 
loss of angular momentum from dynamical friction, it 
will shrink in radius. As a result gas can be transported 
on the dynamical friction timescale to smaller radii. In 
addition, I ffiid that this satellite-disk interaction also 
builds a substantial surface density enhancement at the 
disk edge. This may lead to enhanced star formation at 
the disk edge. This process which I call gas shepherding 
may be a generic feature of gaseous disks around galaxies, 
if a sufficiently massive and dense satellite is available. 

The physics of satellite-disk interactions in galaxies 
or gas shepherding is not just an interesting exercise in 
mathematical physics, but may be important in the fu- 
eling of active galaxies and central starbursts. First, the 
action of dynamical friction on the satellite-disk system 
will shrink radius of the disk, thereby transporting gas 
to smaller radii. This shepherding of gas will continue 
until the gas is forced into the center or the shepherding 
satellite is destroyed. I will discuss this scenario for nu- 
clear fueling and how it compares to competing scenarios 
in this paper. 

I have structured this paper as follows. In fj2l I dis- 
cuss the basic picture of gas shepherding and give a few 
order of magnitude estimates. I estimate the timescale 
for gas shepherding and the expected surface density en- 
hancement. I present the physics of gas shepherding and 
solve numerical models in f}3l I also present an approx- 
imate analytic solution to this problem. I then discuss 
the application of gas shepherding to the feeding of ac- 
tive galactic nuclei and central starbursts in iJH Finally 
I present my conclusions in fJS] 

2. BASIC PICTURE 

The basic picture of gas shepherding is illustrated by 
Figure [T] A gas disk with a radius of rd and surface 
density of S orbits about a spherical mass distribution 
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Fig. 1. — Schematic of gas shepherding. A satelhte with mass, 
Ms, and orbital radius, rg, infalls via dynamical friction and inter- 
acts with a gaseous disk with mass, = nSr^, that sits in a halo 
with total enclosed mass, Monc- 

with enclosed mass of stars or dark matter of Menc(''d)- 
Due to dynamical friction on the stellar background, an 
external satellite with mass Ms slowly spirals in on a cir- 
cular orbit with a radius of rg in the same orbital plane 
of the gas disk. The radial velocity at which the satel- 
lite initially spirals into the gas disk due to dynamical 
friction is Wdf ~ (Ms/MencjVorh, where Vorh is the orbital 
velocity of the satellite. As the satellite approaches the 
disk edge, i.e., rs approaches rd, the satellite excites spi- 
ral density waves at the Lindblad resonances (Goldreich 
& Tremaine 1980), which transfers angular momentum 
from the disk to the satellite. This angular momentum is 
in turn transferred to the background stars or dark mat- 
ter via dynamical friction on the satellite. Since angu- 
lar momentum of the satellite-disk system is continually 
bled, the radius of the satellite-disk system will shrink at 
the shepherding velocity, Vshep- 

I first give some simple order of magnitude estimates of 
this process. The dynamical friction timescale for a satel- 
lite is (Chandrasekhar 1943; Binney & Tremaine 1987) 



^df 



Ms 



dyn, 



(1) 



where t^y^^ ~ fJg ^ is the dynamical time, fJg = 
\J GMonc/rl is the orbital frequency of the satellite, and 
G is Newton's constant. The torque due to dynamical 
friction on the satellite, Tdf, is 



Trif - M^rlVlst-l 



(2) 



while the disk, whose edge sits at rd < r^, torques up the 
satellite (Td). When Tdf ~ Jd, the satellite and the disk 
are well coupled for Ms < Md, i.e., the mass of the satel- 
lite is smaller or comparable to the mass of the disk. 
Setting Td ^ 7d gives 



Md 



1/3 



(4) 



which I define as the "Hill" radius of the disk, rH,d = 
rs(Md/Menc)^'''^-^ This gives the natural scale for the 
radial separation between disk and satellite, i.e., Ts — r^ ^ 

For spiral density waves that damp locally, the radial 
extent of the disk over which the satellite exerts torques 
is also ~ rH,d- Assuming no viscous spreading, the disk 
surface density will be enhanced at the disk edge due to 
the piling up of material from the initial radius of rd,o 
to the current radius of rd. This enhancement due to a 
satellite that sinks a distance ~ rd is 



r-d/rn.d, 

^0 



(5) 



where S is the surface density and Sq is the initial surface 
density. For M^/M^nc ^10^'^ — 10^^, the enhancement 
is a factor of a few to ten. 

3. GAS SHEPHERDING 

3.1. Basic Equations 

To study this problem in greater detail, I begin by 
writing the equations for a viscous disk that is evolving 
under the influence of an external torque. The equation 
of continuity is 



as 1 d{rJ:vr 



dr 



= 0, 



(6) 



where r is the radial coordinate, and Vr is the radial com- 
ponent of the velocity (Frank, King, amd Raine 2002). 
The angular momentum equation is 



djr^n) 
dr 



1 (dT, 



dr 



f 1 ^ 



where VL — ^ GMcncir) jr'^ is the orbital frequency, 
Moncff) is the mass enclosed inside r, and ^visc — 
—^Txr^vYsdSXjdr is the viscous torque, v is the viscos- 
ity (Frank, King, amd Raine 2002). Using equation ([7]), 
I find^ 



The torque due to the excitation of spiral density waves 
in a disk is (Goldreich & Tremaine 1980; Lin & Pa- 
paloizou 1986; Ward & Hourigan 1989; Artymowicz 1993; 
Ward 1997) 



1 



a(r2f]) 



27rrE V dr 
Plugging equation ^ into (O, I find 



f-, (-^visc -^d) 

or 



(8) 



GMt Md 
rs Mpn 



rs - ''d 



(3) 



where Afd = 'I'Sr^ is the mass of the disk, I have pre- 
sumed that radial separation between satellite and disk 
is small compared to the radius, i.e., rs — rd <C rg. 

For now, I assume that there is no internal rearrange- 
ment of angular momentum in the disk, i.e., an inviscid 
disk. Dynamical friction torques down the satellite (Tdf), 



_ 1 d 
dt 27rr dr 



d{r^n) 
dr 



d_ 

dr 



(Tv: 



(9) 



I use quotation marks because disks do not have Hill radii 
in the traditional sense. Rather I have adopted this terminology 
because of the familiar scaling, i.e., the 1/3 power law. 

The angular frequency, C is a function only of the coordinate, 
r, assuming that the mass enclosed, i.e., A/onc is not an explicit 
function of time. Therefore, the first term completely disappears 
from equation (0. 
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For the torque density due to an orbiting sateUite on a 
disk, dT^/dr, I take a form suggested by Ward & Houri- 
gan (1989) (see also Goldreich & Tremaine 1980; Lin & 
Papaloizou 1979a, 1979b, 1986): 

dr 



sgn (r - Ts) 



riSl - a)2 (r - r,)^ 



where /3 is a constant of order unity. For \r 
expand VL — Vl^k, {dVL/dr){r — Ts) to get 

1/3 



(10) 
< r, I 



— « sgn (r - 



(r - rs) 



4 ' 



(11) 



where I have taken a density profile for the enclosed mass 
of the form p (x r~^, so that 

r 



ro 



(12) 



where Monc,o is the enclosed mass at vq. Note that for 
such a mass distribution, the orbital velocity, Vorh = 
yjGMcncfi/rQ, is a constant. Using (fT2|) . I find the fol- 
lowing relations 



dr 



dr 

sgn (r 



= 27r'(;orb^ (Ei/r) 



s)/3- 



GM^Sror^ 



AfencO (j' 

Plugging equations (fTS]) and (fT4|) into I 



5S 



1 (rz^S) 
r dr'^ 



(3 rp GM^ d 
27r Vorhr Monc,o 9r 



), I find 
Er3 



(13) 
(14) 



(r 



(15) 

where I assume r < rg, which fixes the sign. Note that 
the term in equation dH) and ([9]), d{r^il)/dr — Uorb, is a 
constant using the prescribed mass distribution in equa- 
tion (fT2| . I chose the viscosity law suggested by Lin & 
Papaloizou (1986): 

(16) 



The time-rate change of angular momentum of the 
satellite is 

M,^&^ = -T, + T,u (17) 
ot 

where —T^ = — J drdT^/dr. The velocity of infall, i.e., 
the shepherding velocity, is z^shep = rs- Thus, I find 



dr^ 
dt 



forb 



-(3 



Sgn (r 



(r 



cnc,0 



■dr-lnA— ;i8) 



Equations ([TS]) and (|18p constitute a complete set of 
equations which governs the behavior of the satellite-disk 
system. These equations are exactly the same as those 
that govern the migration of protoplanets in protoplan- 
etary nebula (Lin & Papaloizou 1979ab, 1986; Hourigan 
& Ward 1984; Ward & Hourigan 1989; Ward 97; Rafikov 
2002), except with the addition of another term on the 
RHS of equation (fT5|) , which is due to dynamical friction. 

To simplify equations (fTS]) and (HH]), I rescale the vari- 
ables 



t 



r^o i In A 



(19) 
(20) 



where Sq is the initial surface density of the disk, which 
I assume to be constant and VLq is the orbital frequency 
at r = tq. Equations pS]) and (|18p become 



da i^Q 9^ {ra^ 

dt' qflo In Ar dr'^ 
dvs^ ,7T^ f ar^ 
dt' ^ Afenco J {r~r,Y 



d 



(21) 



dr V ^ ''s)'', 
dr - ^ (22) 



where q — Mg/Monc.o and = /3/(27rlnA), and I have 
taken the viscosity to be ^ = vq(j^ (see eq.[16]), where 
= a(/id,o/'"d,o)^''d,o''^orb,^ whcrc Td.o is the initial ra- 
dius, /id,o is the initial vertical scale height of the disk, 
and a dimcnsionlcss viscosity parameter (Frank, King, 
& Rainc 2002). The initial scale height of the disk 

can be written as /id.o/?'d.o = Qoqd/V^,^ where Qq is 
the initial Toomre Q (Toomre 1964) of the disk and 
Qd = Md/Monc,o is the ratio of the disk mass to the 
enclosed mass. Without loss of generality, I now take 
ro = rd^o ~ rsfl to find 

r2 „ ^|2/^^3^ 



dt' q' r dr"^ 

fQ'Qdr^fi d ( ar^ 



r dr \ (r — r, 



drg „ „, f ar^ 



dt 



rdr 



'd.o 



(r - r^Y rs 



(23) 
(24) 



where q' — M^/M^ and a' — aQl /[2h\K). Equations 
(j23p and (|24p are the master equations governing my 
model which I solve numerically in ii3.2l and solve ap- 
proximately in i j3.3l 



3.2. Numerical Results 

Equations ([23|) and ([24|) consists of four dimension- 
less parameters: g'di 9'j oi' and /?'. For the dimension- 
less parameters, /?', which represent the relative power 
between disk tides and dynamical friction, I choose a 
fiducial value for /?' = 1, assuming fiducial values of 
In A ~ 0(1) and (3 ~ 0(1). For a', which represents 
the relative power between viscosity and dynamical fric- 
tion, I take a' — IQ^^ — 10"-'^ to study the effects of 
viscosity. Finally, I set — Md/Monc = 10"'^. 

I solve equations ((23)) and (j24p using standard ex- 
plicit finite-difference methods (Press et al. 1992). For 
good spatial resolution, I select 750 grid points between 
f/''^d,o = to 1.5 with the satellite initially at rg = 2rd.o- 
The satellite begins to fall inward via dynamical friction. 
As the satellite approaches the disk, the disk exerts a 
torque on the satellite, which slows its infall. At the 
same time, the surface density of the disk increases at 
the disk edge. 

I show this evolution for a' = 0.1, (3' ~ 1, and q' — 0.1 
in Figure [51 I start the evolution of the satellite at 

* Using the standard Shakura & Sunyaev (1973) a prescription 
for the viscosity, i.e. u = acsh, where Cs is the sound speed and h 
is the disk scale height, I write Cs = v^rhh/r and I take h = h^ g 
at r = rd,o- 

^ To get this relation, take the Toomre Q parameter for a gaseous 
disk, Q = CsKo/ttCBo (Binney & Tremaine 1987). The epicyclic 
frequency, kq = (4: + 2dlnQo/dlnr)^^'^Qo = y/2Qo. Plugging Cs = 
v^,-t,h/r, I find h/r = {Md/Afcnc)(Q/V2). 
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Fig. 2. — Evolution of the surface density as a function of radius 
for a satellite disk system where g' = 0.1 for a' = 0.1 and = 1. 
The snapshots are separated in At = 0.25ttif q increments. The 
satellite position is marked by the black dot below the surface den- 
sity plot. Each point and each curve correspond to one snapshot. 
The labeled points, A-D, are explained in the text. 

?'s/»'d,o ~ 2. When the satelhte sinks to rs/rd,o — 1-25, I 
set t = and plot the evolution of the satellite-disk sys- 
tem in intervals of At = 0.25tdf,o for a total of 1.75idf,o 
where tdf,o is the dynamical friction timescale at rd,o- At 
t — ^ the disk is undisturbed except for a small amount of 
spreading due to viscosity at the disk edge. The satellite 
position is at r/rd^o ~ 1-25 as noted by the black dot la- 
beled A. I evolve in intervals of At and note that the disk 
satellite system has begun to respond (point labeled B). 
The disk begins to develop a surface density enhancement 
and the satellite slows down a bit. Evolving for another 
At, i.e., the point labeled C, shows a substantial change. 
The surface density increases near the disk edge by a 
maximum factor of « 2. Note that the distance between 
point B and C is substantially smaller than the distance 
between A and B, which results from the backreaction 
of the disk tides. By point D, the density enhancement 
is « 5. Note also that the distance between point D 
and its neighbor is somewhat larger than the other sep- 
arations. This is due to the decrease in the dynamical 
friction timescale compared to the dynamical timescale 
when the mass enclosed decreases relative to the satellite 
mass. 

In Figure [3] and m I show the case for g' ~ 0.01 and 1 
respectively. For = 0.01 (Fig. [3|), the initial infall of 
the satellite is stopped by the disk because of the large 
mass ratio between the satellite and the disk. In addition, 
there is no surface density enhancement at the disk edge. 
For q' = 1 (Fig. |4]), on the other hand, the initial infall 
of the satellite is only slightly slowed because the mass of 
the disk and satellite are the same. In this case, a huge 
surface density enhancement is produced. 

One remarkable aspect of Figures [2] and |4] is that they 
appear to approach quasi-steady state solutions. This is 
especially clear in Figure[21 Motivated by the appearance 
of these quasi-steady state solutions, I find an approxi- 
mate kinematic wave solution in ^ 33. 31 which agrees well 
with these numerical results. 

In Figure [5l I plot the normalized shepherding velocity 
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Fig. 3. — Same as Figure [2] but for g' = 0.01. The larger mass 
of the disk stops the body from infalling toward the center. In this 
case there is no surface density enhancement. 
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Fig. 4. — Same as Figure [2] but for g' = 1. The surface density 
is enhanced over the case where g' = 0.1. In addition the disk has 
a much reduced effect on the infall velocity of the satellite because 
they have equal masses. 

of the satellite disk system, i.e., normalized to the initial 
infall velocity, Wdf.o = "I'd! (?'d,o) (lower plot) and the local 
infall velocity, v^i (upper plot), as a function of for 
= 0.01, 0.1, and 1. I plot it for a' = 0.001, 0.01, and 
0.1. For comparison I plot the local dynamical friction 
infall velocity, fdf (''s)- For > rd,Oi Vshep ~ fdf because 
tidal coupling between the disk and satellite has not yet 
been achieved. Once the coupling has been achieved, i.e. 
''s < '"d.Oi the shepherding velocity depends mainly on 
the mass ratio between the satellite and disk, g', and 
is independent of the viscosity, a' . This highlights the 
fact that dynamical friction dominates the dynamics. I 
also plot an estimate for migration velocity (eq.[lT]) for 
fs < Tdfi ~ ^H,d, i-e., I allow the shepherding satellite to 
fall inside of one Hill radius of the disk to allow a kine- 
matic wave solution to be set up (see ij3.3p . This estimate 
(eq.[4T]) gives remarkably good results, i.e., < 10% for 
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Fig. 5. — Plot of the normalized shepherding velocity as a func- 
tion of Tb for q' = 0.01 (bottom set of curves), 0.1 (middle set), 
and 1 (top set) and a' = 0.001 (short-dashed lines), 0.01 (dotted 
lines), and 0.1 (solid lines) normalized to the initial infall velocity, 
''df,o- The local dynamical infall velocity, v^i (dot-dashed line), is 
also plotted for purposes of comparison. Finally, the estimate for 
''shep f eg. 1411 ) is plotted for rs < q ~ '"H,d to allow a kinematic 
wave solution to be established. 

the fidicial cases of q' — 0.1 and 1. The most remarkable 
aspect of Figure [5] is that while the infall velocity falls 
substantially, a sizeable fraction of the initial infall veloc- 
ity remains. For instance, the infall velocity does not fall 
below 20% of the initial infall velocity for the q' = 0.1 
case. 

The shepherding satellite transports gas from large 
radii to smaller radii by dumping the angular momen- 
tum of the gas into the background stars or dark matter. 
Moreover, this transport is achieved on a dynamical fric- 
tion timescale, which can be fast for sufficiently massive 
satellites. This transport will continue unless the gas is 
consumed via star formation or the satellite is tidally dis- 
rupted by the increasing density of the background (see 



dxs 
dt' 



= Xdfl 



XdflJ J X'^ Xs 



(27) 



where a;s = r's/rH,d and a;d,o = ?'d,o/7'H,d.'' 

I now assume the kinematic wave ansatz: a = a^ix — 
Vshcpt'), where CTw denotes my kinematic wave solution 
for a. Inserting this ansatz into equation (|26p. I find: 



dx q' dx"^ 



d a„ 



ql'^ \Xd,oJ ox x^ 

Integrating once and ignoring the constant of integration, 
I find 



'^shcp 



3a' dal 
2 q' dx ' ^2/3^^ J 



4r.P'{ — ] i> (29) 



where I have cancelled one factor of from each term. 
I now integrate once more to find 



= C- 



I 1/3 

" 9d 



2/3' 
~9~ 



a^d.o 




where C is the constant of integration. The two terms in 
the parenthesis are both negative because x (x r — < 
and Wshop < (inward migration). The first term in 
the parenthesis cuts off the surface density as x becomes 
more negative, i.e., away from the disk edge, while the 
second term cuts off the surface density as it approaches 
the disk edge due to tidal interactions. The peak, i.e., 
where tr^ is maximal, is found from equation ([29]) by 
setting da^jdx = 0. This occurs when* 



X — 3^pcak 



q 

2/3 
.Id 



Xs 
Xdfl 



1 



1/4 



(31) 



Hence, crw(a;poak) = Cmax, so given an estimate for a^t 
I can determine the constant of integration, C, via 



3.3. Approximate Solution 



Motivated by the numerical results of g3T2] I now solve 
equations (|23[) and (|24p using a kinematic wave approach 
inspired by the planetary literature (Hourigan & Ward 
1984; Ward & Hourigan 1989; Ward 1997; Rafikov 2002). 
First I approximate equations (^5)1 and ([M]) to recover 
a simpler set. As the order of magnitude estimate in 
reveals, the typical radial separation for a well coupled 
satellite-disk system is the Hill radius of the disk, rn.d- 
Thus, I now change variables in equations (|23|) and (|24p 
from r to x, where 



(25) 



7'H.dO 



1/3 ' 
Id ^d,0 



1 /3 

where rn.do = 9d ''d.o is the initial Hill radius of the disk. 
Since ru do ^ I take the approximation r — rg ^ r to 
find 



da_ 
'dt' 



'Id 



1/3 



q' dx^ 



q 



a'9d 



1/3 



2 ^shcp-^pcak 



2^ 
~9~ 



Xs 

Xd,o 



2/3 

9d 



peak 



(32) 



In the appendix, I calculate CTmax by balancing the 
torque from the satelite-disk interaction with the torque 
from dynamical friction. I find the following cubic equa- 
tion (eq. [M3] V 



2/3' a„ 



peak 



Xd,Q 



2 / 

a 9d 3 



a^d,o 

Xa 



= 0, (33) 



^ To aid in the derivation of equation l|26|l and l|27|l . it is helpful 
to consider which terms will change rapidly for small x and which 
terms will change slowly. For instance, the term (r — rs)~* — > a;"* 
changes rapidly with small x, whereas r — > r = rs -|- rjj d^; ~ ''s 
changes very slowly. Hence, I made the following series of substi- 
tutions r — > rs and r — rg — > rjj d^;- 

* To clarify the signs in equation JSTJ, note that 2: < 0, i.e., 
the satellite is at larger radii than the gas and v^hcp < 0> the 
satellite-disk system moves inward. 
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which I solve numerically. This gives an estimate for 
Cmax and hence, C, given Wshop- Therefore, to complete 
the calculation, I now give an estimate for V ghcp ■ Before 
I do so, I note that the RHS of equation pO|) can be 
negative. To join this kinematic wave solution, CTw, onto 
the background solution, I set 



if cr^ > and x > Xpcak 
if cr^ > 1 and x < Xpoak 

if cr^ < and X > Xpeak 
if cr^ < 1 and X < Xpeak 



(34) 



I now proceed to calculate Ushop- I estimate the torque 
using equation ([2]). The satellite raises tides on the disk, 
which do not operate over the entire disk, i.e. the tides 
are most acute at the edge of the disk due to the (r — 
rs)~^ scaling of the torque density. Rather they act on 
the accumulated mass at the disk edge, Mcdge, which I 
estimate to be: 



edge 



1 



^d,cdgc 
fd,0 



(35) 



where r^ cdgc is the edge of the disk.^ To estimate rd.cdg 
for equation (|35p . I presume the following fit 



^■d.cd 



do, 



(36) 



where ip is a constant of order unity and is determined 
below. The total angular momentum of the satellite-disk 
system over which these torques are effective is 



L = Uorb {Msrs + Modgc^'d) 



(37) 



The difference between rd and Tg is of order the Hill ra- 
dius of the disk, rn^d = rd{Md/Mcnc)^^^ (see As the 
satellite falls inward and builds a large surface density 
enhancement, I expect that the Hill radius of the disk 
will change, i.e., Sq ^ S « Mcdgc/47rrdrH,d, where I as- 
sume that the mass of the edge is spread over a radial 
thickness of rn.d- Therefore, I define an effective Hill 
radius of the disk by 



r-H,eff _ / 7r£:(rd)r ^ ^ 

rd V Mcnc 

_ m,d / 

ro \2TTMdr}i,d 



edge 



1/3 



(38) 
(39) 



The disk edge is separated from the satellite by ^ rH.ofr- 
Thus, Td « (1 — f7?'H.cff/''d)''s, where 77 is a constant of 
order unity. Equating the time derivative of L fea.[37j) 
and (121), I find 



GM2 



L — V, 



orb 



Ms + M, 



edge 



1 



V- 



? ^H,eff 



where I ignore the time derivatives on Mcdgc- 
Ms Ms 



Vshc 



' Voih 



f^shepj 

(40) 
This gives 

(41) 



Me„c(r,) Ms + Modge (1 - ?7m,off/rd) 



Despite it crudeness, this velocity estimator gives reason- 
ably good agreement (< 10%) when I select ry = 1 and 

^ In the nomenclature of this paper, ^dgo should really be r^. 
However, I will introduce this variable as I will fit for this and 
^'H.cff • 




Fig. 6. — The surface density enhancement S/Sq as a function of 
(r — ''s)/rH,d for the approximate solution (dashed hnes) and the 
numerical solution (solid lines) for a satellite disk system where 
q' = 0.1 for a' = 0.1 and /3' = 1. The snapshots are separated 
in At = 0.5tdf increments and the corresponding position of the 
shepherding satellite (from top to bottom) is rs/r^ Q = 0.95, 0.75, 
and 0.65. 

= 2 via trial and error. I plot the comparison in Figure 

El 

Equations and dH]) constitute a 

complete analytic solution once the evolution of the disk 
can be modeled as kinematic wave, i.e., after the satellite 
sinks down a few rjj.d- This approximate solution agrees 
well with the exact numerical solution as shown in Figure 
[Hfor q' = 0.1. The agreement for q' — I and q' — 0.01 
is not as good. For q' = 1, the assumption that Td « 
Tdf breaks down and the agreement can be improved by 
estimating Td « Td{/{1 + q'). For q' — 0.01, kinematic 
wave solutions are not established because the infall of 
the satellite is stopped. 

Finally, from the estimate for t^shep given in equa- 
tion (|4T|) and Figure El I note Wshep is a factor of a 
few larger than (Ms/Md)wdf, where Vdf is taken at rd,o- 
Hence I estimate the timescale for shepherding to be 
ishcp f{Md/Ms)tdi or 



^shcp 

where / < 1 and tdf ^ 



MdM^nc tdyn 

Ms Ms In A' 

(Me„c/M,)tdyn/lnA. 



(42) 



4. THE FUELING OF CENTRAL STARBURSTS AND 
ACTIVE GALACTIC NUCLEI 

As the satellite-disk system should continually shrink 
due to dynamical friction, it raises the interesting ques- 
tion of whether or not such a mechanism may be deliver 
gas into the nuclear regions of galaxies. The fueling of 
central starbursts and active galactic nuclei remains an 
open question. It is generally believed that the fueling of 
central starbursts and active galactic nuclei (AGNs) re- 
quires the delivery of gas from large radii to small radii 
(for a review see Shlosman, Begelman, & Frank 1990; 
hereafter SBF90). Gas transport is mediated at small 
radii (~ 1 pc) by a thin accretion disk. At much larger 
radii (> 0.1 — 1 kpc), galactic bars can efficiently trans- 
port gas. 
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How is gas transported in the gap between 0.1 — 1 
kpc to ^ 1 pc? Galactic bars, while efficient at large 
scales, stop at the inner Lindblad resonance (ILR), which 
is at ~ 0.1 — 1 kpc.^" Thin accretion disks, which are 
sufficient at small scales, would become gravitationally 
unstable and would fragment to form stars outside of 
~ 1 pc (Paczynski 1978; Kolykhalov & Sunyaev 1980; 
Shlosman & Begelman 1987; Goodman 2003; Thompson, 
Quataert & Murray 2005). 

There are several possible modes by which this trans- 
port may be accomplished. First the gas is transported 
by large-scale gaseous gravitational instabilities, i.e., the 
"bars within bars" model (Shlosman, Frank, & Begel- 
man 1989). Secondly, the gas may form a thick tur- 
bulent marginally stable accretion disk powered by star 
formation (Paczynski 1978; Kolykhalov & Sunyaev 1980; 
Shlosman & Begelman 1987; Goodman 2003) via radia- 
tion pressure (Goodman & Tan 2004; Thompson et al. 
2005) or supernova (Wada & Norman 2002). Third, if 
the gas exist in the form of molecular clouds, collisions 
between them may lead to episodic feeding of the nu- 
clear regions, i.e., the chaotic accretion scenario (Krolik 
& Begelman 1988; SBF90; Hopkins & Hernquist 2006; 
King & Pringle 2007; Nayakshin and King 2007). Fi- 
nally, the fueling of these central regions by be one of 
brute force, where major or minor mergers drive gas from 
the disk of the galaxy to the nuclear region (Hernquist 
1989; Bekki & Noguchi 1994; Mihos & Hernquist 1994; 
Hernquist & Mihos 1995) 

In the "bars within bars" scenario, Shlosman, Frank, 
& Begelman (1989) argued that if a gravitationally un- 
stable gas disk is sufficiently massive such that it fulfills 
the Ostriker-Peebles criterion (Ostriker & Peebles 1973), 
it forms a secondary gaseous bar which funnels gas in- 
ward on a dynamical time. Numerical simulations con- 
firmed the viability of "bars within bars" (see for instance 
Friedh & Martinet 1993; Levine et al 2007). It demands 
a substantial mass in gas, Mgas ~ Afenc, where M^nc is 
the enclosed mass. While such large gas concentrations 
are observed (see for instance Jogee et al. 2005), de- 
livered into the nuclear region by large scale bars, it is 
unclear if this mechanism drives the bulk of AGN fuel- 
ing and quasar activity. Namely, the observed correlation 
between large-scale bars, which transports gas into the 
central region of galaxies, and Seyfert activity is statisti- 
cally marginal (see the review by Ho 2008 and references 
therein; see the review by Combes 2003 and references 
therein; but also see Laine et al 2002) . In addition, there 
is no evidence that Seyferts have a higher fraction of nu- 
clear bars (see Knapen 2004 and references therein). 

On the other hand, if the disk is gravitationally un- 
stable but lacks sufficient mass to fulfill the demands of 
"bars within bars" , the gas will fragment and form stars. 
The process of star formation may feed back on the disk 
to provide vertical support, maintaining marginal sta- 
bility, i.e., (5^1, where Q is the Toomre Q (see for 
instance Shlosman & Begelman 1989; Goodman 2003; 

^'^ As pointed out by the referee, ILRs may or may not exist, 
depending on the nuclear mass concentration. However, these bars 
still become inefficient at small scales, i.e., typically 10% of the size 
of the bar (Shlosman, Frank, & Begelman 1989). 

As pointed out by the first referee, however, the nuclear bars 
may be very short lived and therefore, may be very difficult to 
observed. 



Goodman & Tan 2004; Thompson et al. 2005). These 
thick marginally stable disks may have a sufficiently large 
scaleheight such that accretion via viscous processes is 
possible before all the gas is consumed. 

Third, perhaps the gas exist primary in the form of 
molecular clouds. Collisions between these molecular 
clouds provides something akin to a viscosity (Krolik & 
Begelman 1988; SBF90). These clouds would then feed 
the nuclear region in a chaotic fashion (Hopkins & Hern- 
quist 2006; King & Pringle 2007; Nayakshin and King 
2007). Typically, the viscosity due to cloud-cloud col- 
lision appears to be insufficient for fueling the central 
regions from hundreds of parsecs. However, it may be 
necessary around the sphere of influence of the central 
black hole where large scale instabilities may not oper- 
ate (SBF90; see fgH) 

Fourth, the fueling of central nuclear activity need not 
be due to gas in the ISM, but gas liberated from dis- 
rupted stars. Hills (1975) proposed that the tidal disrup- 
tion of stars (Rees 1988) may be a means by which active 
galaxies may be fueled. These rates of tidal disruption 
are enhanced in mergers of galaxies (Roos 1981), which 
links the three phenomenon of merging galaxies, active 
galactic nuclei, and tidal disruption of stars. Alterna- 
tively, a pre-existing gas disk may capture stars plunging 
through it, subsequently disrupting these stars, which re- 
plenishes the disk. This gas may then viscously accrete, 
fueling the central engine (Miralda-Escude & KoUmeier 

2005) . 

Finally, mergers between galaxies may drive gas to- 
ward the nuclear regions of galaxies where they fuel nu- 
clear activity (Hernquist 1989; Bekki & Noguchi 1994; 
Bekki 1995). This possibility may fuel ultraluminous in- 
frared galaxies and quasar activity as modeled, for in- 
stance, by Springel et al. (2005) and Hopkins et al. 
(2005) and has significant observational support (see for 
instance Sanders et al 1988). However, active galactic 
nuclei (AGN) activity appears to be due to more quies- 
cent phenomenon (Li et al 2006; Hopkins & Hernquist 

2006) . For instance, Li et al. (2006) explored the clus- 
tering properties of AGN in Data Release 4 of the Sloan 
Digital Sky Survey (SDSS) and concluded that interac- 
tions are unlikely to be the main culprit in the fueling 
of AGNs. Minor mergers are another possibility (see for 
instance Mihos & Hernquist 1994; Hernquist & Mihos 
1995; Tanguchi & Wada 1996; Taniguchi 1997, 1999). 

4.1. Fueling via Gas Shepherding 

Gas shepherding is another scenario by which gas is 
tranported from hundreds of parsecs. The gas shepherd- 
ing scenario, which I present is as follows: large scale 
bars drive gas toward the central regions of a galaxy 
where they collect in a nuclear ring, i.e. the gas tran- 
sitions from Xi to roughly circular X2 orbits (Binney & 
Tremaine 1987) at a radii of a few hundred parsecs to 1-2 
kiloparsecs (Buta & Combes 1996; Knapen 2004). This 
gaseous ring becomes gravitationally unstable and forms 
giant molecular clouds which in turn collapse to form 
star clusters in coplanar orbits (see for instance, Maoz et 
al. 2001; Jogee et al. 2002). If the star cluster formed is 
sufficiently massive, it may shepherd in the remnant gas 
in the gaseous ring from a few hundred parsecs into the 
center. 

I first compare the timescale of shepherding and vis- 
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cous processes. For an a viscosity, the viscous tiniescale 
is ivisc ~ {r/h)^tdyn (Frank, King, & Raine 2002). 
For a marginally stable disk, i.e., Q ^ 1, h/r ^ Md/Mcnc- 
Therefore the ratio between the shepherding timescale, 
^shop: and the viscous timescale, ^visci is 

Uep fa Md fMdV 

2 



10 



-3 



0.01 



//lnA\ fMd/M, 



0.1 



(43) 



These timescales become comparable to each other for 
thicker disks, i.e. h/r ~ M^/Mcnc ~ 1, or lower mass 
satellites, i.e., Md/Mg ~ 100 or a combination of the 
two. Thus when shepherding operates, it dominates over 
viscous processes for low mass disk and large satellites. 
Therefore, for this scenario to work, the star cluster that 
is formed must be sufficiently massive compared to the 
mass of the gas disk. I would expect it to be at least 10% 
or greater in order for shepherding to be efficient based 
on equation (HS)) and on my calculations in SjH] 

If gas shepherding operates, what ultimately hap- 
pends to this shepherded gas? First, if this gas is in a 
marginally stable disk, it may fragment into stars at the 
disk edge (see for instance Thompson et al. 2005; Levin 
2007; Chang et al. 2007; Nayakshin, Cuadra, & Springel 
2007; Van der Ven & Chang in preparation). Secondly, 
this gas may be shepherded into a radius where the shep- 
herding satellite is tidally disrupted. I now discuss these 
two possibilities below. 

4.2. Disk Fragmentation and Star Formation 

Marginally stable disks with sufficiently rapid cooling 
(icooi < 3<dyn) will tend to fragment (Gammie 2001). In 
the galactic context this fragmentation will lead to star 
formation, which may feed back on the disk, maintain- 
ing marginal stability. This feedback may be due to mo- 
mentum driving (Murray, Quataert, & Thompson 2005; 
Thompson et al. 2005), cosmic ray feedback (Socrates, 
Davis, & Ramirez- Ruiz 2008), or supernovae driven tur- 
bulence (Wada & Norman 2002). A surface density en- 
hancement of a few to ten would enhance local star for- 
mation. I take the standard Kennicutt law (Kennicutt 
1998) for the star formation rate to be 



(44) 



where is the star formation rate per unit area, E 
is the gas surface density, 77 ~ 0.01 is the empirically 
determined efiiciency. Equation (|44|) suggests that the 
star formation rate would also be enhanced by the same 
factor as that of the surface density. 

Equation (j44|l implies that the timescale for the gas to 
fragment into stars is 

J] 

^ ^ ^ rr^tdyn- (45) 

The ratio between the timescale for shepherding and the 
timescale for the gas to fragment into stars is 



t 



shcp 



J_ 

In A 



Md 



(46) 



1 



0.01 



/ .//In A 



0.1 



Md/M, 



Menc/Md 

100 



(47) 



Thus, sufficiently large satellites can shepherd a consider- 
able fraction of the gas into the central regions before the 
gas completely fragments into stars. For typical numbers 
(Menc ^ 10^ M0), this implies a satellite mass ~ 10^ M©. 

For less massive satellites, the disk would completely 
fragment into stars. Thus, the evolution would likely 
be different. As the satellite continues to fall inward, it 
captures these newly-formed stars into mean-motion and 
secular resonances, which increases the star's eccentric- 
ity and inclination. Indeed, Yu, Lu, & Lin (2007) have 
suggested that this exact process may be responsible for 
the distribution of the young stars in the Galactic Cen- 
ter (GC). In this case, an intermediate-mass black hole 
(IMBH) or dark cluster falls into a disk of stars in the 
same orbital plane. The stellar disk is heated up as the 
stars are captured into the previously-mentioned reso- 
nance. Some of these stars are forced to migrate with 
the infalling IMBH, while others are forced into close 
encounters. Some of these stars that are subjected to 
these close encounters may be ejected out of the GC 
at large velocities, producing hypervelocity stars (HVSs) 
or hypervelocity binaries (HVBs) (Lu, Yu, & Lin 2007). 
Hence stars, whose semi-major axis is larger than the 
current radius of the infalling IMBH would be excited 
into a torus-like structure, while stars interior to the ra- 
dius of the infalling IMBH would remain in a cold disk. 
Such a dynamical distribution may result from the effects 
of shepherding, but at much larger scales, i.e., ~ 100 pc 
rather than ~ 0.4 pc as in the case in the GC. 

4.3. Destruction of the Shepherd 

If this gas does not fragment into stars (which is the 
case for sufficiently massive satellites), gas shepherding 
may transport a substantial amount of gas from ~ 0.1 
kpc to smaller radii. As long as the satellite remains 
coherent, there is no limit to how far gas may be shep- 
herded. This requirement that the satellite be sufficiently 
massive (~ 10^ M©) limits the number of possible sys- 
tems that can serve as a shepherd. 

A super star cluster (SSC) is one such candidate. SSCs 
range up to mass of ~ 10^ — 10^ Mq with a radius of par- 
sees (see for instance McCrady & Graham 2007; Hagele 
et al. 2007) so they are sufficiently massive to trans- 
port gas before it completely fragments into stars. How- 
ever, as they sink deeper into the bulge, tidal forces may 
disrupt them, stopping the shepherding process. The 
SSC is tidally disrupted when the enclosed background 
density penc oc r~^ is equal to their mean stellar den- 
sity, pssc ~ -Afgsc/aisC' where agsc is the radial size 
of the SSC. Plugging in a few example numbers, I find 
that pssc - lO6(Afssc/lO6M0)(assc/lpc)-3M0pc-3. 
For typical numbers, i.e., Monc ~ 10^ Mq at 100 pc, 
Pcnc ~ 10^ (r/l pc)~^ Mq pc~'^ so the radius where the 
satellite is disrupted is 

^"'^^^H sxIO^Mq J (^j 
which would also be where the gas disk is shepherded 
down to. Note that I have ignored the effect of the 
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SMBH, which begins to become important at such radii. 
The densities and density structures of these SSCs are 
not very well constrained. McCrady and Graham (2007) 
presents a very detailed study of the size and mass distru- 
bution of super star clusters in the nearby starbursting 
galaxy M82. The mean half power radius of these SSCs 
is 1.8 pc (McCrady and Graham 2007, their Figure 7), 
and there are quite a few clusters around 5 x 10^ M0. 
By equation this gives a tidal disruption radius of 
about 10 pc. Lets consider two examples in more de- 
tail. One example is SSC L in M82 which has a virial 
mass of « 4 X 10^ Mq, a projected half-light radius of 
« 1.5 pc, and a tidal disruption radius of rtid ~ 3pc. 
Another example is SSC 11, which has a virial mass of 
« 4 X 10^ M0, a projected half light radius of w 1.1 pc, 
and a tidal disruption radius of sa 6 pc. In another study, 
Hagele et al (2007) observed five star forming regions of 
the circumnuclear ring of NGC 3351 and found star clus- 
ter masses between 1.8 — 8.9 and radii of 1.7 — 4.9 pc. 
This would bracket tidal disruption radii between a few 
parsecs to 25 parsecs. Hubble space telescope studies by 
Maoz et al (1996) suggest that the star clusters formed 
in circumnuclear rings have sizes < 5 pc. Hence, SSCs 
appears have sufficient mass and density to shepherd gas 
into ten to tens of parsecs. 

The destruction of a SSC may not be the result of tidal 
processes, but instead to internal processes, i.e., stellar 
mass loss. Thus, SSCs may have a short lifetime in anal- 
ogy with stellar clusters in our galaxy, i.e., open clusters 
and OB associations. For typical galactic clusters, they 
are disrupted because gas ejection from the formation of 
HH regions and supernova (Lada and Lada 2003). How- 
ever, SSCs are thought to be the progenitors of globular 
clusters, whose lifetime is ^ 10^° years, and thus their 
lifetime may be substantially longer. The vulnerability of 
SSCs to these internal processes strongly depend on their 
stellar initial mass functions, i.e., they must have a suffi- 
cient number of low mass stars (Meurer et al 1995). Top- 
heavy IMFs are vulnerable to the disruptive processes of 
stellar mass loss and other processes (Chernoff & Wein- 
berg 1990). Young, dense stellar clusters, the closest 
local analogues to SSCs, show a significant population of 
low mass stars. For instance, the 2 x lO"' Mq (Walborn 
et al 2002) stellar cluster, R136 in the Large Magellanic 
Cloud has a large population of low mass stars down to 
0.1 M0 (Brandl et al 1996). In addition, McCrady et al. 
(2003) found SSC 9 was consistent with a normal IMF, 
though SSC 11 appeared to have a top heavy IMF. 

Aside from SSCs, the nucleated remains of an accreted 
satellite galaxy may also be the shepherds that herds that 
gas into the central regions of the galaxy for sufficiently 
small mass ratios. For mass ratio between the mass of 
the satellite galaxy's remnant and the dynamical mass 
is of order 10% or greater, long range tidal interactions 
may disrupt a central disk before the shepherding sce- 
nario outlined aboved can be set up and would be akin 
to the minor merger scenario studied previously (Hern- 
quist 1989; Bekki & Noguchi 1994; Mihos & Hernquist 
1994; Hernquist & Mihos 1995; Tanguchi & Wada 1996; 
Taniguchi 1997, 1999). Using the nucleated remains of 
satellite galaxies as shepherds is interesting, but a de- 
tailed study of this scenario is left to future work. 

4.4. The Final Ten (or Tens of) Parsecs 



As I argue in SSCs can shepherd gas to within 

10 to 10s of parsecs of the center. However, due to tidal 
disruption of the shepherds, gas shepherding is less ef- 
fective at these small radii as a mechanism to drive gas 
further inward. A similar problem exists for the "bars 
within bars" scenario. Once gaseous bars funnel the gas 
within the gravitational sphere of influence of a black 
hole, it becomes increasingly difficult to have a sufficient 
amount of mass to form nested bars (SBF90 their Figure 
2). As a result, another means of gas transport between 
ten to tens of parsecs must operate for the gas shep- 
herding and "bars within bars" scenarios, i.e., the "10 
parsec" problem. One possibility, outlined by SBF90, 
are cloud-cloud collisions. At the deposition radius, the 
gas becomes marginally self-gravitating and fragments 
into molecular clouds whose collisions with one another 
serves as a source of viscosity which drives material into 
the central regions. Typically, the viscosity scales like: 

( C for C « 1 

«ci-ci = I ^-1 ^ 1 , (49) 

where C is the filling fraction of the clouds (see SBF90). 
It is unclear if this is more effective than magnetic (see 
the review by Balbus 2003) or gravitationally driven vis- 
cosities (Gammie 2001). 

The effects of cloud-cloud collisions on the velocity dis- 
persion of the clouds is not known. The internal struc- 
ture, i.e., magnetic fields, will play a huge role in de- 
termining the efficiency of cloud-cloud coUisional driven 
viscosity (see for instance Krolik & Begelman 1988). An- 
other way by which this final "10 pc problem" can be 
solved is via a disk which is supported either by ra- 
diation pressure (Thompson et al. 2005), supernova 
driven turbulence (Wada & Norman 2002; Kawakatu & 
Wada 2008), or cosmic ray pressure (Socrates, Davis, & 
Ramirez-Ruiz 2008). The gas disk or ring which results 
from the deposition of gas at ten to tens of parsecs may 
form a marginally stable disk which is supported by star 
formation/supernova which can then viscously accrete. 
One-zone calculations by Thompson et al. (2005) suggest 
that a star formation rate of ~ 1 Mq yr~^ is sufficient to 
support a marginally stable disk at ten to tens of par- 
secs (see Thompson et al. 2005, their Figure 5),^^ while 
also allows for a nontrivial accretion rate onto the central 
black hole. I note that the strengths of the "bars within 
bars" or gas shepherding scenarios and that of the star- 
burst supported disk scenario are complementary to one 
another. Inside of tens of parsecs, modest star formation 
rates can support a radiation pressure supported disk 
rather than requiring 200 — 6OOM0yr^^ star formation 
rates from 200 pes inward (again see Thompson et al. 
2005, their Figure 5). Recent observation of NGC 3227 
by Davies et al. (2006) support this notion of star for- 
mation in a gas disk at small radii. However, the current 
star formation rate of ~ 0.1 M0yr~^ appears to be in- 
sufficient to support the vertical scale of the disk. Other 
forms of pressure support, i.e., cosmic ray pressure (see 
Socrates, Davis, & Ramirez- Ruiz 2007 their appendix E), 
may be important in this context. 

5. DISCUSSION AND CONCLUSIONS 

I note, however, that they presume some sort of global torque 
via spiral arms or the like in their calculations. 
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I have calculated the action of an infalling satellite on a 
coplanar gaseous disk. The satellite falls inward because 
of dynamical friction and couples to the gaseous disk via 
disk tides. Due to the unrelenting loss of angular momen- 
tum, the satellite-disk system shrinks in radius. Thus, 
gas is transported from larger radii to smaller radii. I 
calculate the structure of the evolving disk as a function 
of radius, both numerically and analytically and show a 
significant density enhancement at the disk edge. This 
density enhancement is due to the transport of gas ini- 
tially at larger radii to a ring that is rn.d thick radially. 

The fate of this shepherded gas is not known. I pre- 
sented two possible outcomes: 1. the gas will fragment 
and form stars or 2. the gas is shepherded into the 
nuclear region. In the latter case, I suggest that the 
gas shepherding scenario, which I have presented may 
be one means by which the "100 parsec" problem may 
be resolved. For this to occur, gas must be shepherded 
sufficiently rapidly such that the gas is not completely 
consumed in star formation, which requires a sufficiently 
massive satellite (~ 10^ M©). I argue that such massive 
satellites exist on the form of SSCs. 

In the model presented, the physics of star forma- 
tion is a disk environment has not been fully addressed. 
Observationally, galactic star-forming disks appear to 
be marginally unstable to fragmentation, i.e., Q ^ 1 
(Thompson et al 2005) and their star formation rate fol- 
lows the Kennicutt law (Kennicutt 1998). These two 
points suggests that some sort of feedback is responsible 
for maintaining this careful balance between fragmen- 
tation and heating. How this feedback operates is not 
well understood. As discussed previously, various models 
have invoked radiation pressure (Thompson et al. 2005), 
supernova feedback (Wada & Norman 2002), and cosmic 
ray pressure (Socrates et al. 2008) from star formation. 
Alternatively, the turbulent pressure support may re- 
sult from gravitational instability, i.e., gravitoturbulence 
(Gammie 2001). It is unclear that in the presence of gas 
shepherding if this equilibrium can be maintained at the 
edge of the disk where the surface density is enhanced. 
If it does not, this edge is likely to fragment completely 
into stars on a few dynamical times. It is also unclear 
if the Kennicutt law would continue to hold at the disk 
edge. Additional insights from future research on self- 
gravitating star-forming galactic disks will be needed to 
help formulate a more self-consistent calculation of the 
gas shepherding scenario, which will be needed to fully 
address the fate of the shepherded gas. 

Assuming that the gas is not fully consumed in star for- 
mation, it would be shepherd down to 10 to 10s of parsecs 
(see M.4p . At that radii cloud-cloud collisions (Krolik & 
Begelman 1988; SBF90) or a starburst supported disk 
(Thompson et al. 2005) would allow the material to be 
viscously accreted. 

The initial satellite inclination and/or eccentricity is 
likely to significantly affect gas shepherding. In this pa- 
per, I have studied a special case where the orbit of the 
satellite is in the plane of the disk and approximately cir- 
cular for simplicity. An eccentric satellite would have its 
eccentricity damped via the excitation of spiral density 
waves. If the eccentricity is below rn^^/r^, i.e., the satel- 
lite stays within one "Hill radius" of the disk of its outer 
edge, then the shepherding calculation presented should 
be a reasonable estimate of the shepherding timescales 



and velocities. 

Compared to a coplanar satellite, an inclined satel- 
lite would couple less well to the gaseous disk. For a 
sufficiently large inclination, the gas disk and satellite 
would be effectively decoupled. However, its inclination 
could be sufficiently rapidly damped via induced bend- 
ing waves in the gas disk in addition to spiral density 
waves (see for instance Ostriker 1994; Terquem 1998). 
Once its inclination is damped below the critical value, 
where the disk and satellite become well coupled, the 
shepherding scenario presented in this paper would likely 
follow. I would expect this critical inclination to be of 
order sini = ru.d/^'d, i-c, the satellite stays within of 
order a "Hill radius" of the disk above and below it. 

Extending the present calculation to satellites of mod- 
est inclination would be fruitful and would allow a study 
of using the nucleated cores from minor mergers as shep- 
herding satellites (see ^4.3p . As discussed in Sj4l niinor 
mergers may drive nuclear activity. The initial inclina- 
tion of minor mergers is arbitary, but their inclination is 
damped as they interact with the large scale stellar and 
gaseous disk. As a result, they may approach the shep- 
herding scenario presented in this paper. Such a study 
is a subject for future work. 

In addition, dynamical friction would also act on the 
satellite-induced spiral density waves. Dynamical fric- 
tion of the extended mass perturbation, i.e., the satellite- 
induced spiral density waves, may enhance the transfer 
of angular momentum between the gas disk and back- 
ground stars (Tremaine & Weinberg 1984). A study of 
these effects, while fruitful, is beyond the scope of this 
work. 

The physics of gas shepherding is also important in 
other areas as well. For instance, it may shape radial 
distribution of star clusters in nuclear gas rings. In some 
of these nuclear ring systems, the location of the star 
clusters is exterior to the gas ring from which they are 
presumably formed. For instance, NGC 4314 clearly 
show the star clusters exterior to the gas (Benedict et 
al. 2002). Martini et al. (2003) found that the eight 
galaxies with strong bars all have star formation exterior 
to the dust ring in their sample of 123 galaxies. How gas 
shepherding drives this morphology is the subject of a 
forthcoming work. (Van der Ven & Chang, in prepara- 
tion). 

I thank N. Murray for his early encouragement and 
useful discussions. I thank A. Socrates for his advice 
and continuing encouragement. I thank R. Levine and 
A. Kratsov for making a early copy of their manuscript 
available. I thank E. Quataert for useful discussions and 
a careful reading of this manuscript. I also thank G. 
Van der Ven and E. Chiang for useful discussions. I 
thank M. Jones for a careful reading of this manuscript. I 
thank the anonymous referees for very useful comments 
which improved the presentation of this paper. I also 
thank the Kavli Institute for Theoretical Physics and 
the Max-Planck Institute for Astrophysics for their hos- 
pitality during the initial and final stages of this project, 
respectively. I gratefully acknowledge the support of the 
Miller Institute for Basic Research. This research was 
supported in part by the National Science Foundation 
under Grant No. PHY05-51164 



Gas Shepherding 



11 



APPENDIX 
CALCULATION OF ctmax 

I now discuss the calculation ol CTmax and hence derive the constant ol integration C . 1 first plug equation (|30p into 
equation ([27| to find 



\ 



c 



q' 



2/3' f X, 



I 1/3 

"9d 



( \\\dx 2:^,0 



2/3 

9d 



(Al) 



I make the approximation that the torque due to the satelhte-disk interaction and the torque from dynamical friction 
nearly cancel each other out, i.e., x^^dxsjdt' <^ 1. Hence I find. 



2/3' 



C 



q' 



1/3 



2/3' f Xs \ q' f 1 \\ dx f Xd.o 

2/3 

9d 



It is helpful to break the integral in equation (|A2p into two parts. I rewrite the integral, /, as 

C^pcak 



J = 



+ 



(A2) 



(A3) 



where Xmax < a^s is defined as where ay„{x > Xmax) = 0, i.e., the disk is cutoff by disk tides above this radius, and 
consider each integral in turn. The first integral, Ii , is 



C 



\ "'9d 



1/3 



2(3' ( X, 



9 Vxdn 



2/3 
^d 



1 



X'^ I X' 



dx 



(A4) 



At 

^peak; — ^max and it declines as x becomes more negative, i.e., away from the satellite, because of the Vshep^ 
term in the square root. However, the torque density declines even faster, so I may approximate Cw = fmax for the 
purposes of this integral, i.e.. 



-dx 



1 CTn 



The second integral, I2, is 



l2 = 



C 



2(3' f X, 



\ /I 



/3 



3 9 V^d,0 



2/3 
^d 



1 \ \ da; 



x-^ I x^ 



(A5) 



(A6) 



Here, the disk torque increases as x approaches aJmax, but CTw declines due to the x~^ term in the square root. As the 
linear term Vshc-pX is overwhelmed by x~^ as x approaches smaller absolute values i.e. x a;max, I ignore the linear 
term and approximate this integral by 



'C- 



g'2 2/3' / X., 



a'qd 9 \Xd 



1 \dx 



(A7) 



Changing variables to 



and performing the integral, I find 



X. 



y^C- 



q'^ 2/3' / X, 



a'qd 9 V^d^o 

' ^'2 / ^ 2 



dy- 



1 



2^^ 
3 a'qd \Xd,oJ x 



rdx. 



C- 



2(3' ( X, 



I 1/3 
"9d 



3 9 \xdfi 



2/3 

9d 



3/2' 



^pcak 



(A8) 
(A9) 

(AlO) 



At Xpoak the term in brackets, which I identify as my approximation for CTw, is crj^iax; while at Xmax, the term in brackets 
is 0. Thus 1 find 

\ Xs J q' 



/2 max 



(All) 
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Chang 



/^/.+I.4^+/3-(^)'^aL.. (A12) 



Hence equation (|A3|) becomes 



Plugging this resuh into equation (|A2p gives a cubic equation for an 
which I may solve numerically. 
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